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Abstract 

We determine theoretically the relation between the total number of protons Np and the mass 
number A (the charge to mass ratio) of nuclei and neutron cores with the model recently proposed 
by Rufhni et al. (2007) and we compare it with other A'p versus A relations: the empirical 
one, related to the Periodic Table, and the semi-empirical relation, obtained by minimizing the 
Weizsacker mass formula. We find that there is a very good agreement between all the relations 
for values of A typical of nuclei, with differences of the order of per cent. Our relation and 
the semi-empirical one are in agreement up to A ^ 10*; for higher values, we find that the two 
relations differ. We interprete the different behaviour of our theoretical relation as a result of the 
penetration of electrons (initially confined in an external shell) inside the core, that becomes more 
and more important by increasing A; these effects are not taken into account in the semi-empirical 
mass-formula. 

1 Introduction 

It is well known that stable nuclei are located, in the Nn-Np plane (where Nn and Np are the total 
number of neutrons and protons respectively), in a region that, for small values of A'p, is almost a line 
well described by the relation Nn = Np. 

In the past, several efforts have been made to explain theoretically this property, for example with the 
liquid drop model of atoms, that is based on two properties common to all nuclei: their mass densities 
and their binding energies for nucleons are almost indipendent from the mass number A = Nn + Np 
[3]. This model takes into account the strong nuclear force and the Coulombian repulsion between 
protons and explains different properties of nuclei, for example the relation between Np and A (the 
charge to mass ratio) . 

In this work we derive theoretically the charge to mass ratio of nuclei and extend it to neutron cores 
(characterized by higher values of A) with the model of Ruffini et al. [2 , [3] . We consider systems 
composed of degenerate neutrons, protons and electrons and we use the relativistic Thomas-Fermi 
equation and the equation of /3-equilibrium to determine the number density and the total number of 
these particles, from which we obtain the relation between Np and A. 



2 The theoretical model 

Following the work of Rufhni et al. [2], [5], we describe nuclei and neutron cores as spherically 
symmetric systems composed of degenerate protons, electrons and neutrons and impose the condition 
of global charge neutrality. 
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We assume that the proton's number density np{r) is constant inside the core (r < Rq) and vanishes 
outside the core (r > Rc)- 

Mr)-(Sr)(^iRc~r), (1) 



where Np is the total number of protons and Rc is the corc-radius, parametrized as 

Rc = A— (2) 

We choose A in order to have p ^ pN, where p and pjq are the mass density of the system and the 
nuclear density respectively (pjv — 2.314 ■ g / cm~^) . 
The electron number density ne(r) is given by 

^e(r) = ^[pf(r)]', (3) 

where (r) is the electron Fermi momentum. It can be calculated from the condition of equilibrium 
of Fermi degenerate electrons, that implies the null value of their Fermi energy (r) 

ef (r) = y^[pf(r)c]2 +m2c4 - rrieC^ + Vc{r) 0, (4) 

where Vc{r) is the Coulomb potential energy of electrons. 
From this condition we obtain 



P. 

hence the electron number density is 



^{r) = - v/l42(^)_2mec2K(r), (5) 



The Coulomb potential energy of electrons, necessary to derive ne{r), can be determined as follows. 
Based on the Gauss law, 14 (r) obeys the following Poisson equation: 

VV,(r) = -47re2K(r) - np(r)], (7) 

with the boundary conditions Vc{oo) ~ 0, Vc(0) = finite. Introducing the dimensionless function 
defined by the relation 

Vc{r) = -hc^, (8) 
r 

and the new variable x = 7'&^^ = r (^ip;^^ , from cq. ([7]) we obtain the relativistic Thomas- Fermi 
equation 



1 d\{x) J 1 4 ^..2...^ ^ ...m3/2- 

- " Td[Xc ~ X) 



(9) 



3x dx^ [ A3 ' ' Qtt 

The boundary conditions for the function xi^) 

X(0) ^ 0, x(oo) = 0, (10) 

as well as the continuity of xi^) its first derivative x i^) s-t the boundary of the core. 
The number density of neutrons n„ (r) is 

^n{r) = -^K{rf, (11) 
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where Pni^) is the neutron Fermi momentum. It can be calculated with the condition of equilibrium 
between the processes 

+ p ^ n + Ve] (12) 

n ^ p + er + i7e^ (13) 

Assuming that neutrinos escape from the core as soon as they are produced, this condition (condition 
of /^-equilibrium) is 

efM+<M=e^M- (14) 

Eq. (fT4|) can be explicitly written as 



\J[Pp{r)cY + m^d^ - TTipC^ - Vc{r) = \f [Pni''")^]'^ + m^C^ - rUnC^ ■ (15) 

3 Np versus A relation 

Using the previous equations, we derive rie(r), n„(r) and np{r) and, by integrating these, we obtain 
Ne, Nn and Np. We also derive a theoretical relation between Np and A and we compare it with the 
data of the Periodic Table and with the semi-empirical relation 

that, in the limit of low A, gives the well known relation Np = A/2 [J. 

Eq. ()16|) can be obtained by minimizing the semi-empirical mass formula, that was first formulated 
by Weizsacker in 1935 and is based on empirical measurements and on theory (the liquid drop model 
of atoms). 

The liquid drop model approximates the nucleus as a sphere composed of protons and neutrons (and 
not electrons) and takes into account the Coulombian repulsion between protons and the strong nuclear 
force. Another important characteristic of this model is that it is based on the property that the mass 
densities of nuclei are approximately the same, indipendcntly from A [T]. In fact, from scattering 
experiments it was found the following expression for the nuclear radius Rn'- 

RN=roA^/^ (17) 
with rg = 1.2 fm. Using eq. (|17|) the nuclear density can be write as follows: 

ArriN SArriN irriN 



where rriN is the nucleon mass. From eq. (|18|) it is clear that nuclear density is indipendent from A, 
so it is constant for all nuclei. 

The property of constant density for all nuclei is a common point with our model: in fact, we choose A 
in order to have the same mass density for every value of A; in particular we consider the case p ~ pat, 
as previously said. 

In table Q are listed some values of A obtained with our model and the semi-empirical mass 
formula, as well as the data of the Periodic Table; in fig. ^ and ^ it is shown the comparison 
between the various Np — A relations. 

It is clear that there is a good agreement between all the relations for values of A typical of nuclei, 
with differences of the order of per cent. Our relation and the semi-empirical one are in agreement up 
to A 10^; for higher values, we find that the two relations differ. We interprete these differences as 
due to the effects of penetration of electrons inside the core [see fig. ([3])]: in our model we consider a 
system composed of degenerate protons, neutrons and electrons. For the smallest values of A, all the 
electrons are in a shell outside the core; by increasing A, they progressively penetrate into the core [5]. 
These effects, which need the relativistic approach introduced in [2j, are not taken into account in the 
semi-empirical mass formula. 

We also note that the charge to mass ratio become constant for A greater that 10^; in particular, 
it is well approximated by the relation Np = 0.026 A [see fig. (|4])]. 
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Figure 1: The Np — A relation obtained with our model and with the semi-empirical mass formula, the Np 
relation and the data of the Periodic Table; relations are plotted for values of A from to 200. 



= A/2 




Figure 2: The Np — A relation obtained with our model and with the semi-empirical mass formula and the Np = A/2 
relation; relations are plotted for values of A from to 10*. It is clear how the semi-empirical relation and the one 
obtained with our model are in good agreement up to values of A of the order of 10*; for greater values of A the two 
relation differ because our model takes into account the penetration of electrons inside the core, which is not considered 
in the semi-empirical mass formula. 




Figure 3: The electron number in units of the total proton number Np as function of the radial distance in units of 
the core radius Rc, for different values of A. It is clear that, by increasing the value of A, the penetration of electrons 
inside the core increases. Figure from Ruffini et al. [2]. 



Figure 4: The Np — A relation obtained with our model and the asymptotic limit Np = 0.02671 





Aa/ 


ApT 


AsE 


5 


10.40 


10.811 


10.36 


10 


21.59 


20.183 


21.15 


15 


32.58 


30.9738 


32.28 


20 


44.24 


40.08 


43.72 


25 


56.17 


54.938 


55.45 


30 


68.43 


65.37 


67.46 


50 


120.40 


118.69 


118.05 


70 


176.78 


173.04 


172.54 


90 


237.41 


232.038 


230.79 


110 


302.18 


271 


292.75 


150 


443.98 




427.73 


200 


644.03 




617.56 


250 


869.32 




831.63 


300 


1119.71 




1071.08 


350 


1395.12 




1337.23 


450 


2019.48 




1955.57 


500 


2367.77 




2310.96 


550 


2739.60 




2699.45 


600 


3134.28 




3122.83 


10^ 


6.9-10^ 




8-10^ 


10* 


2.0-10'^ 




3.45-10« 


10^ 


3.0-10« 




3.38-10'J 


10*^ 


3.4-10^ 




3.37-10^^ 


10'^ 


3.7-10« 




3.37-10^'^ 


10^" 


3.9-10^^ 




3.37-10^'* 



Table 1: DifFerent values of Np (column 1) and corresponding values of A from our model {Am, column 
2), the Periodic Table {Apt, column 3) and the semi-empirical mass formula {Ase, column 4). 
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4 Conclusions 



In this work we have derived theoretically a relation between the total number of protons Np and the 
mass number A for nuclei and neutron cores with the model recently proposed by Ruffini et al. [5], 

We have considered spherically symmetric systems composed of degenerate electrons, protons and 
neutrons having global charge neutrality and the same mass densities (p Pn)- By integrating the 
relativistic Thomas-Fermi equation and using the equation of /3-cquilibrium, we have determined the 
total number of protons, electrons and neutrons in the system and hence a theoretical relation between 
Np and A. 

We have compared this relation with the empirical data of the Periodic Table and with the semi- 
empirical relation, obtained by minimizing the Weizsacker mass formula by considering systems with 
the same mass densities. We have shown that there's a good agreement between all the relations 
for values of A typical of nuclei, with differences of the order of per cent. Our relation and the 
semi-empirical one are in agreement up to A ~ 10^; for higher values, we find that the two relations 
differ. We interprete the different behaviour of our theoretical relation as a result of the penetration 
of electrons (initially confined in an external shell) inside the core [see fig.(l3|)], that becomes more and 
more important by increasing A; these effects, which need the relativistic approach introduced in [2], 
are not taken into account in the semi-empirical mass-formula. 
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